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We study lepton asymmetry evolution in plasma of the early Universe before the electroweak 
phase transition (EWPT) accounting for chirality flip processes via Higgs decays (inverse decays) 
entering equilibrium at temperatures below Trl — 10 TeV, Tew <T < Trl. We solve appropriate 
kinetic equations for leptons and Higgs bosons taking into account the lepton number violation due 
to Abehan anomalies for right and left electrons and neutrinos in the self-consistent hypercharge 
field obeying Maxwell equations modified by the contribution of the Standard Model of electroweak 
interactions. The violation of left lepton numbers and the corresponding violation of the baryon 
number due to sphaleron processes in symmetric phase is taken into account as well. Assuming the 
Chern-Simons wave configuration of the seed hypercharge field, we get the estimates of baryon and 
lepton asymmetries evolved from the primordial right electron asymmetry existing alone as partial 
asymmetry at T > Trl. One finds a strong dependence of the asymmetries on the Chern-Simons 
wave number. We predict a nonzero chiral asymmetry = /Xep^ — /Xel 7^ in this scenario evolved 
down to the EWPT moment that can be used as an initial value for the Maxwellian field evolution 
after EWPT. 
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I. INTRODUCTION 

The nature of the initial fields that seed subsequent 
dynamo for observed galactic magnetic fields is largely 
unknown [T|,ll. It might be that seed fields are produced 
during the epoch of the galaxy formation, or ejected 
by first supernovae or active galactic nuclei. Alterna- 
tively to this astrophysical scenario the seed fields might 
originate from much earlier epochs of the Universe ex- 
pansion, down to the cosmological infiation epoch Q. 
There are first observational indications of the pres- 
ence of cosmological magnetic fields (CMF) in the inter- 
galactic medium which may survive even till the present 
epoch 

It is well known that Maxwellian CMF might arise dur- 
ing the electroweak phase transition (EWPT) from mass- 
less (long-range) hypercharge fields Y"^ existing in pri- 
mordial plasma before EWPT Q. Note that long-range 
non-Abelian magnetic fields (corresponding to, e.g., the 
color SU(3) or weak SU(2) groups) cannot exist because 
at high temperatures the non-Abelian interactions induce 
a "magnetic" mass gap ~ g'^T. 

The Faraday equation which governs evolution of mag- 
netic (hypermagnetic) fields in the early Universe de- 
pends crucially on the helicity parameter a which is 



* maxim.dvornikov@u sp.br] 
t ,semikoz@yandex.ru| 



a scalar owing to the parity violation in the Standard 
Model (SM) and having different forms before and after 
EWPT. We remind that the standard magnetohydrody- 
namic (MHD) parameter, which is generated by vortices 
in plasma ctMHO ~ (v- (V x v)) is pseudoscalar according 
to the parity conservation in QED plasma. Obviously in 
the isotropic early Universe such vortices are absent, at 
least at large scales we consider here. 

Before EWPT in the symmetric phase of primordial 
plasma the helicity parameter ay results from the Chern- 
Simons (CS) anomaly term in the SM Lagrangian 0] 
Lcs — (3'^A^eR/47r^)(BY-Y), which violates parity. Here, 
g' = e/ cos^w is the Uy{1) gauge coupling, By = V x Y 
is the hypermagnetic field, and /igR is the chemical poten- 
tial for right-handed electrons (positrons). The polariza- 
tion origin of the CS term Lcs was elucidated in Ref. Q 
as an effect of comoving right electrons and positrons 
along By which have opposite spin projections on By 
and populate the main Landau level with slightly differ- 
ent densities due to fien, 7^ 0. 

While such CS term vanishes in broken phase at T < 
Tew 0, there appears a similar polarization effect 
due to weak interaction of neutrinos (antineutrinos) with 
polarized electrons (positrons) given by the axial vector 

(A) (A) 

force Fcr = — VV^ arising from the parity violating 
part of weak interactions between neutrinos of the flavor 
a = e, fj,,T and charged leptons, 

yi-4)(x,t)-GFci^)(M(-).5j(''")(x,i)). 

This axial vector force acts only on those polarized elec- 
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trons and positrons which contribute to the partial mag- 
netization M*^°'^ = //B sign((T) uq^'B/B with the corre- 
sponding densities noa of electrons (tr = — ) and positrons 
(cr ~ +) at the main Landau level. Here Gp is the 

(A) 

Fermi constant, Cq = =f0.5 is the axial coupling in the 
Wcinberg-Salam model, the upper (lower) sign stays for 
electron (muon, tau) neutrinos, and /ib is the Bohr mag- 
neton. Note that the neutrino current density asymmetry 
^j(i>a) — — jp„ is the polar vector and the magneti- 
zation M*^"^"^ is the axial vector. Thus, the weak axial- 
vector force separates electric charges causing the rel- 
ative drift velocity and the corresponding electric cur- 
rent that leads to the generation of an additional elec- 
tromagnetic field component resulting from weak inter- 
actions E„oak ^ ckB Gp. Unfortunately the a- 
helicity parameter produced by this mechanism is very 
small. Moreover, it depends on such additional pa- 
rameter as a neutrino gas inhomogeneity scale A^^''jjj, 

a ~ Gf{T / x''^^.^){Sni, /n^) [l^]. It should be mentioned 
that the new opportunity for the Maxwellian CMF gen- 
eration found in recent works [ll|, [T2| seems to be very 
intriguing. 

In Ref. [ll| it was shown that the evolution of 
Maxwellian magnetic fields in a primordial plasma at 
temperatures T > 10 MeV is strongly affected by the 
quantum chiral anomaly proportional to the difference 
of right-handed and left-handed electron chemical po- 
tentials A/i(i) — ^Lf,^ — /ieL evolving in a self-consistent 
way. Such a difference defines the new helicity parame- 
ter anow(^) — acmA/i(t)/7r(Tcond in the modified Faraday 
equation that governs evolution of Maxwellian fields just 
after EWPT at temperatures 10 MeV < T < Tew (see 
comments below in Appendix IX)) . Here acm = 1/137 is 
the fine-structure constant and Ucond is the plasma con- 
ductivity. 

The goal of the present work is a more careful analy- 
sis of the initial chiral anomaly parameter A ^(Tew) aris- 
ing in hypermagnetic fields before EWPT and accounting 
for the Higgs (inverse) decays which change chiralities of 
right- and left-handed electrons (positrons). While the 
sphaleron processes are always switched off in the Higgs 
phase, in the symmetric phase at T > Tpw, conversely, 
we should take them into account since lepton and baryon 
numbers are violated even in the absence of hypermag- 
netic fields. The plan of our paper is the following. In 
Sec. [H] we formulate equilibrium conditions in the sym- 
metric phase of primordial plasma and comment on the 
corresponding set of chemical potentials. In the main 
Sec, mil we derive kinetic equations for appropriate asym- 
metries of right- and left-handed electrons (positrons) in 
the presence of Abelian anomalies for fermions, as well as 
the kinetic equation for the Higgs boson asymmetry ac- 
counting for both inverse and direct Higgs decays. Such 
an extension of the kinetic approach, considered earlier 
in our paper [l^, allows us to correct the value Afi{TE,w) 
accounting for the Higgs boson asymmetry evolution. In 
Sec. nil Al we check conservation laws and calculate the 



baryon asymmetry of the Universe (BAU) through lep- 
togenesis in hypermagnetic fields. In Sec. |IV] we solve 
kinetic equations for all lepton and Higgs asymmetries 
both analytically, neglecting hypercharge fields, and nu- 
merically, including hypermagnetic fields. In Sec. |V] we 
analyze the evolution of the chiral anomaly parameter 
Afi{t) in the symmetric phase down to t ~ tpw- Our 
results are discussed in Sec. I VII 

In Appendixl^we interpret and compare two quantum 
mechanisms producing a-helicity parameter for magnetic 
and hypermagnetic fields: the chiral anomaly leading to 
the ancw-helicity parameter in the Faraday equation [ll| , 
and the Chern-Simons anomaly for hypercharge fields in- 
terpreted as a polarization effect in plasma caused by the 
hypermagnetic field itself Q . In Appendix |B] we derive 
kinetic equations for the lepton and Higgs asymmetries 
used in the main Sec. IIIII In Appendix [C] we give some 
formulas for the lepton number violation due to 't Hooft's 
anomaly in non- Abelian fields in order to explain how the 
sphaleron processes influence the left lepton kinetics. 

II. EQUILIBRIUM IN THE SYMMETRIC 
PHASE OF ELECTROWEAK PLASMA AND 
THE CHIRAL ANOMALY PROBLEM 

The question of how large the chiral anomaly param- 
eter Afi{t) = fien it) — MclC*) could bc before EWPT is 
important as an input for the generation of Maxwellian 
magnetic fields after EWPT. In the SM plasma consist- 
ing of quarks, leptons, and one Higgs doublet ip"^ — 
{^p'-^\(p'^^'>), with the chemical potential being in Bose 
distribution /iq = ^,^(0) — one can expect the 

chemical equilibrium in symmetric phase given by the 
relation 

A^ea - Mel = "MO- (2.1) 

Equation (|2.ip corresponds to Higgs decays and inverse 
decays in reactions eLCR O (p'^"-' and v^eB. ^ (p^'^\ Here 
for the SM doublet = {v^ei,) we use the equality of 
chemical potentials in Fermi distributions /igL = Mi^^ ■ 

For the case of the global equilibrium in the absence of 
hypermagnetic fields, similar reactions with Higgs bosons 
obey analogous relations both for other lepton genera- 
tions, jiij^ — jjLi^ = —Ho, I — ^j.,T, and for quarks, e.g., 
Mml - fJ-dn = Mo, Mmr - Ai«L = Mo, in reactions qu^ Qd^ ^ 
and (7«a9"L ^ f''^^ correspondingly (T^. ITa|. 

We consider below only one generation with the low- 
est Yukawa coupling of fermions with Higgs bosons hg = 
^/2me/v = 2.94 x 10~^. Thus right electrons enter the 
equilibrium with left particles through Higgs (inverse) 
decays in the expanding Universe in the last instance. 
This is because of the high rate of chirality flip reac- 
tions Frl ^ h^T, which becomes faster than the Hubble 
expansion H ^ , Frl > H, at temperatures below 
Trl ^ 10 TeV. This fact is important in scenarios where 
the generation of BAU proceeds through the leptogenesis, 
and a primordial BAU is stored in right electrons that 
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are being protected from washing out by sphalerons all 
the way down to Trl- Suggesting such scenario, the au- 
thors of Ref. [iB] supposed that such value Trl is close 
to the temperature at which the sphaleron effects fall 
out of the equilibrium, and therefore it is possible that 
the CR may not be transformed into eL soon enough for 
the sphalerons to turn them into antiquarks, and thereby 
wipe out the remaining BAU. 

In this scenario the global equilibrium [13, [IE] fails, and 
five (=5) remaining chemical potentials describe equilib- 
rium in a hot plasma before EWPT: three fii for the 
three global charges B/3 — Li = const, where i = 1,2,3 
enumerates generations in SM, fiy for the conserved hy- 
percharge (global {Y) = 0), and /ie^ for right electrons br, 
with the conservation of their lepton number dfj_jl^^ = 
unless T > Trl .6] • Then, if one assumes the presence of 
large-scale hypercharge fields in the symmetric phase, 
which are progenitors of Maxwellian fields in the broken 
phase, the number of right electrons is not conserved be- 
cause of the Abelian anomaly [l3| 

where and are, respectively, the C/y(1) hyper- 
charge field strengths and their duals, and Fr = —2 is 
the hypercharge of the right electron. 

There are no asymmetries of left leptons and Higgs 
bosons in this scenario, /Iel = Mo 0, and the chiral 
asymmetry (|2.1|) reduces to — ^e^- For such sce- 
nario with a nonzero cr asymmetry alone sphaleron 
washing out BAU is absent all the way down to EWPT. 

In a broadened scenario with nonzero left lepton asym- 



metries = ^u'^ 7^ 0, where Ci = f^a/T, appropri- 
ate for the stage T < Trl [1, [Ts], we somehow vio- 
late the equilibrium described in Ref. [^ by five chemi- 
cal potentials for five globally conserved charges. Nev- 
ertheless, it can lead only to an additional factor of 
the order one ca ^ 1 that describes the dependence of 
riL — {rie^ — rig^) = ^el^^/S ^ on five global charges in 
primordial plasma. For instance, rewriting the canonical 
Abelian anomaly for the left doublet Lj — [v]^ , cl), 

«. = -^v^'"^ ^L = -i, (2.3) 

in the form dCeJdt = -CA(6g'Vl67r2T3)(EY • By), we 
put below ca = 1 simplifying the solution of our kinetic 
equations for the lepton and Higgs boson asymmetries. 
Note that assuming a nonzero left particle asymmetry 
7^ 0, we should take into account the sphaleron pro- 
cesses violating lepton and baryon numbers. The compe- 
tition of such processes with hypermagnetic field contri- 
bution through Abelian anomaly is one of the interesting 
questions touched upon in the present work. 

III. KINETICS OF LEPTONS AND HIGGS 
BOSONS IN HYPERMAGNETIC FIELDS 

In Ref. we forced the presence of zero Higgs asym- 
metry n^(o) — n^(o) = = 0, ^0 = considering 
leptogenesis with the inverse decays only, crCl — > <^^°\ 
eRJ^e — >■ t/''-"-' , etc. Now let us consider both inverse Higgs 
decays and direct Higgs decays. The system of kinetic 
equations for leptons accounting for Abelian anoma- 
lies (|2.2p . (|2.3I) . and sphaleron processes for left leptons 
takes the form 



-^(Ey . By) + 2rR, - L... - llV^^^ZV^ 



(-) 



dt A^-^s"-^ ^' ' i 2s 

for decays (inverse decays) crcl O ip^^^ and eRj/g ^ tp 



dt IDTT^S 2 2s 

for brCl ^ ^^^^ , as well as 

^ ^ -(Ey-By) ^L^L -t- Frl |ie^ — i 



for env^ ^ ip'-+\ (3.1) 



Here Lb — (rib — nb)/s is the lepton number, h — 
eRjeLi^^e: s ^ 27r^5*T^/45 is the entropy density, and 
g* = 106.75 is the number of relativistic degrees of 
freedom. The factor of 2 in front of the rate Frl in 
the first line takes into account the equivalent reaction 
branches. We also included Higgs decays with the rate 
Fd Frl/2. The probability Fsph = Ca^T is given by 



sphaleron transitions decreasing the left lepton numbers 
and therefore washing out BAU, where aw — 5^/47r ~ 
1/137 sin^6'w = 3.17 x 10~^ is given by the gauge cou- 
pling g — e/ sin^w in SM and is the Weinberg angle. 
The constant C ~ 25 is estimated through lattice cal- 
culations (see some comments on 't Hooft's anomaly in 
Appendix [Cl and Chap. 11 in Ref. [111). Of course, for 



4 



the left doublet = 



This system is completed by the kinetic equation for 
the Higgs bosons independent of Abelian anomaly inher- 
ent in fermions 19'| 



d 
di 



(3.2) 



Note that the rate of Higgs decays (inverse decays) 
coincides with the rate of a lepton pair production (an- 
nihilation) having opposite sign since the creation of a 
pair is followed by the disappearance of a Higgs boson 
and vice versa. 

In kinetic Eqs. (13. ip and (13.21) we used the rate of all 
inverse processes [16], which is twice bigger than for the 
decay ones Trl = 2rD, 



rH,= 5.3 X 10-3;.? f^'lV^ 



1 



\T J 



2t 



EW 



(3.3) 



This rate vanishes just at the EWPT time x = 1, 
where the variable x — t/t-^w = (2ew/T)^ is given 
by the Friedmann law. Here = 2.94 x 10~^ is the 
Yukawa coupling for electrons, Fq = 121, and ra'^{T) = 
2DT^{1 — T^y^/T^) is the temperature dependent effec- 
tive Higgs mass at zero momentum and zero Higgs vac- 
uum expectation value. The coefficient 2D sa 0.377 for 
mg(T) is given by the known masses of gauge bosons 
niz and mw, the top quark mass m*, and a still prob- 
lematic zero-temperature Higgs mass, which is estimated 
as TTiH ~ 125 GeV (see Ref. [2(^). Of course, the chi- 
rality flipping rate exists after EWPT. However, that 
rate is due to electromagnetic processes at T < Tew, 
Fern — CKcm("^e/3^^)T' whcu particlcs (electrons and 
positrons) acquire the nonzero mass mg. 

The detailed derivation of kinetic Eqs. p.ip and p.2|) 
accounting for chirality flip processes (without Abelian 
anomaly and sphaleron transitions) is given in Ap- 
pendix |B] 

Let us rewrite Eqs. p.ip and p.2[) using the asym- 
metries LgR = Cej^T^/ds, = CeL^^/Ss and (n^(o) - 
n^m)/s = foT^/Ss as 



Ey • By + 2Frl( — Cen + Cel ^ Co), 



dt 27r2T3 

+ rRL(feR - Cel + Co), 

33'^ T, r^ph 

HiY • riy ?eL 



d^^L 



dt Stt^T^ 



dCo 
dt 



+ rRL(CeR - Cel + Co) 
=rRL(-CeR + Cel - Co)- 



(3.4) 



The third equation for neutrinos is excess since — ^ej^ . 
Thus, we have three equations for three chemical poten- 
tials instead of the two in Ref. UM ■ Note that we should 



have d^o/dt < for our initial conditions Cer(^o) > 
and CelI^o) = Co(^o) = resulting in the negative chem- 
ical potential for the boson doublet (p^ = {(p^~^\ (p^'^''), 
/io < 0, as it should be. 

Below we simplify the Abelian anomaly contribution 
~ (Ey • By) considering, as in Ref. [l3|, the simplest 
configuration of hypermagnetic field: CS wave = 
y(t) sin fco^, Yy — Y{t) cos koz, = Yq = 0. Using 
the generalized Ohm's law [8| 

Ey = —V X By + ?7yV x By — ayBy, 

where 7]y = ((Xcond)-^ is the magnetic (hypermagnetic) 
diffusion coefficient, ay is the hypermagnetic helicity pa- 
rameter arising due to the polarization of electroweak 
plasma [sl il3l]. 



ay 



5'^(MeR 



MeJ2) 



47r2(7cond 



f^cond 



lOOT, 



(3.5) 



and we get the pseudoscalar (Ey • By) entering the 
Abelian anomaly as 



(Ey • By) =w(V x By) • By - q^yB^ 



100 



fco 

T 



47r2 



Cr 



2 



(3.6) 



Here we substituted (V x By) - By = koB^{t) for the CS 
wave, where By{t) = koY{t) is the hypermagnetic field 
amplitude. 

Using the notations Unix) — 10^^eR(a;), y^ix) — 
10^Cel(2;), and yo{x) — 10''^o(2;), as well as accounting 
for Eq. p.6p . the system (|3.4I) can be rewritten in the 
form that is analogous to Eq. (3.4) in Ref. [l3| (without 
contribution of neutrinos, which is identical to that of 
left electrons) 



da; 



dyL 
dec 



dyo 
da; 

Here 



To ^(2/R 



B 



(0) 



1020 G 



Box'^'-Ao{yn + f 



5.6 X lO^C 



_ro(l-a;) 
2Vx 



yR - yo) 




2,3/2gV(a:^) 



Bo = 25.6 



IO-^Tew 



Ao = 77.6, 



(3.7) 



(3.8) 



are constants chosen for hypermagnetic fields normalized 
on 1020 G. 
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The function e"^^^^ is given by the hypermagnetic field 
squared 



By{x) 



H 2 



(0) 

Y J 



(3.9) 



We also substituted the hypermagnetic field Byit) = 
k[)Y{t) found as the solution of the modified Faraday 
equation [HI, ^ for the CS wave ^ 



Byit) =B^^ exp I [aYit')ko - k^T^y it')]dt' 
= i?fexp(3.5^ 



0.1 



IO-^Tew 
ivK + yL/2) 



fco 



IO-^Tew 



Ax' 



(3.10) 



Note that we do not consider here a negative value of the 
wave number k^ < Q that is allowed as well and could 
lead to the lepton number violation via Abelian anomaly 
proportional to the pseudoscalar p.6p . (Ey • By) ^ 
k^Y'^^t) < 0. This is because the case fco < corresponds 
to the decay of the hypermagnetic field p.lOp instead of 
a real instability evolving in MHD plasma for fcp > 0. 

We choose initial conditions at xq = 10"'* or at Tq — 
Trl when Higgs (inverse) decay becomes faster than the 
Hubble expansion Frl > H, 



ynixo) = 10 ^, Vhixo) = yo(a;o) = 0. 



(3.11) 



Such conditions correspond to the right electron asymme- 
try (,eii{xa) = 10^^° chosen at the level of baryon asym- 
metry. 



A. Conservation laws and BAU in hypermagnetic 
fields 

One can see from kinetic Eq. (|3.1I) that in the absence 
of hypercharge fields the total lepton number is not con- 
served due to sphaleron transitions washing out the left 
lepton number ALe/At = L^^ -\- L^^ + L^l = — TsphieL- 
The baryogenesis arises through the leptogenesis due to 
the conservation law B/3 — — const, where B — 
(riB — n^)/s. Accounting for Abelian anomalies in sys- 
tem p.ip . such baryogenesis is possible, B ^ 0, since the 
hypermagnetic fields raise the lepton number and BAU 
as well die/dt|sY#o > 0, AB/dt\BY^o > 0. This growth 
proceeds opposite to the competing s pha leron influence 
erasing B (compare in Ref. [13| where we ne- 

glected sphaleron transitions) . 

Three global charges are conserved {Si = const). 



B 



3 f"^^ 



<^2, Lr 

' 3 



(3.12) 



as well as I/gR = <^r well above Trl, T ^ Trl. If the ini- 
tial BAU differs from zero, B{to) ^ 0, and if we assume 
the absence of lepton asymmetries for the second and 
third generations all the way down to TeWi = Lr — 0, 
we find that the relation 62 — S3 — B{xo)/3 is valid only 
for the initial time. From the first conservation law in 
Eq. p.l2p one finds the change of BAU B{t) at temper- 
atures T < Trl- This change obeys the relations 



B{t) 



Leit) 



B(to) 



Len (to) 



Si. 



If, for simplicity, we assume the zero initial BAU B{to) — 
or 52,3 — 0, then finally we get the conservation law 

B{t)/3'^Le{t)=-LeM- 

Thus, in the present scenario, BAU sits in hypercharge 
fields that are decreasing due to the sphaleron processes, 
as follows from the sum of kinetic Eq. (|3.1I) : 



B{t) =3 r 



AL,^,{t') , dLeJi') , dL.^Ji') 



At' 

I2\ l-t 



At' 



At' 



At' 



|,j/(E..B.,r. 



A^ 

s 



3 / TsphLei^At' . 

Jto 



(3.13) 



Using the first equation in the system p.7p . where the 
hypermagnetic term comes from the Abelian anomaly 
~ (Ey • By), one obtains from Eq. ()3.13|) the baryon 
asymmetry in the following form: 



B{x) =2.14 X IQ-^ I Ax'i 
(1-x') 



Ax' 



'■^ Ax' 
128C / ^yL(x'). 



(3.14) 



The baryon asymmetry p.l4p for different values of the 
parameter Bq ~ 25.6 x (fco/10~^TEw) or for different CS 
wave numbers fco is shown in Fig. [T] Notice that for very 
small fco ^ fcmax — IO^^Tew, the role of the hypermag- 
netic field, which feeds BAU growth, becomes negligible 
since By ~ fco- As a result, sphaleron transitions wash 
out BAU, or they diminish the Abelian anomaly lepto- 
genesis effect in such a way that BAU can be even neg- 
ative at the EWPT time B{tEw) < [see curve B{t) in 
Fig. [T](b) plotted for the parameter So = 2 x lO^'^ that 
corresponds to fco = 7.8 x 10~^fcmax]- 



IV. CHEMICAL EQUILIBRIUM WITH AND 
WITHOUT HYPERMAGNETIC FIELDS 

Neglecting the hypermagnetic field contribution and 
sphaleron transitions and using the initial condi- 
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(a) (b) 

FIG. 1. The baryon asymmetry B{t) versus t/tEw for By '' ~ 10^^ G. (a) The baryon asymmetry for Bo = 2.1 x 10^^ (sohd 
hne) and Bo ~ 25.6 (dashed hne). (b) The baryon asymmetry for Bo = 2 x 10"'^. 



tion p. lip , we easily find the solutions of kinetic equa- 



tions ((XT)) 












yhix) = 


4 ^ 




yo{x) = 


4 ^ 




$(a;) = 


-4ro 





Note that yo < 0, as it should be for a boson chemical 
potential. Obviously the chemical equilibrium (|2.ip is 
settled soon due to huge negative $ ~ — 4ro = —484, 

2/R-yL+yo = zjR(xo)e*(^) ^0, (4.2) 

that happens somewhere at a; > Xcq — 10^^ at the 
temperature T = Tew/V^ - ^ ^eV before EWPT, 
T > Tew (see right panel of Fig. 

The numerical solution of the system p.7p accounting 
for sphalerons and in the presence of By 7^ for the 
particular case of CS wave configuration is shown on the 
left panel in the same Fig. [51 We see that opposite to 
the case in Eq. (14. 2p . in the presence of hypermagnetic 
fields and accounting for sphaleron transitions, the chem- 
ical equilibrium between leptons and Higgs's (|2.ip never 
exists: the sum ?/r — j/l + yo even grows in the symmet- 
ric phase when t — > iEw- The shorter CS wavelength 
(e.g., for Bq = 25.6 if we use the maximum wave number 
fco = 10~^Tew), the larger values (j/r — J/l + yo) evolve 

up to ^EW- 

All curves in Fig. [5] start from the same initial condi- 
tion ynixo) - Uhixo) +yo{xo) = ynixo) = lO"*^ that cor- 
responds to the initial right electron asymmetry = 
10~i° close to the BAU value we expect at the EWPT 
time X = 1. Thus, the violation of lepton numbers in 



external fields through Abelian anomalies and sphaleron 
transitions leads to a violation of the chemical equilib- 
rium (j2.ip existing in primordial plasma when perturba- 
tive reactions (decays, scattering, etc.) are taken into 
account only. 

V. CHIRAL ANOMALY PARAMETER 
iVR - yL) / IN ELECTRO WEAK PLASMA 
BEFORE EWPT 

The temporal evolution of the chiral anomaly param- 
eter — j/L = 10^(A/Li/T) is shown in Fig. |3l There is 
a strong dependence on the scale A = k^^ for the cho- 
sen configuration of the hypermagnetic field: the shorter 
CS wavelength, the bigger the chiral anomaly parame- 
ter will be. For the maximum acceptable wave number 
ko = fcmax = 10~^Tew (^0 = 25.6), the chiral anomaly 
parameter A/x/T is close to 5 x 10~^ assumed in Fig. 1 
of Rcf. 11] as a maximum initial value of the chiral 
anomaly parameter just after EWPT. However for longer 
CS wavelengths [see dashed line in Fig. ^h) plotted for 
fco = 7.8 X 10~^fcmax], such an initial value will occur 
at a negligible level (A/i/T ~ 10~^) which can crucially 
change the results of Ref . [ll| . Note that for the strongest 
CS field amplitude By = fcmax ^^(Oi the case of the 
maximum Abelian anomaly leptogenesis effect, the left 
lepton asymmetry t/l grows from the beginning due to 
Higgs (inverse) decays and then changes sign, t/l < 0, 
which is allowed for fermions, cf.. Fig. [3l^c). 

Let us explain qualitatively the growth of the chiral 
anomaly shown in Figs.[3l^a) and[31[b). One can simplify 
the kinetic equations for f and in the system p.4p 
decoupling them. For this purpose we neglect the Higgs 
boson asymmetry = 0- We also omit the asymmetry 
of left leptons = in the first line of Eq. p.4p . and 
the right electron in the second line of Eq. p.4p ^ej, = 0. 
For example, from the first equation in the system (13.41) . 
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(a) 



(b) 



FIG. 2. The function 2/r — i/l + yo versus t/t-Ew for lepton asymmetries ya ~ 10*(/ia/r)- (a) The numerical solution of the 
system ([33 for B-^"' = 10^^ G. The solid line corresponds to Bo = 2.1 x 10"^ and the dashed line to Bo = 25.6. (b) The 



analytic expression for j/r — yL + yo given by Eq. (14.2 



substituting the pseudoscalar value (Ey • By) for the CS 
wave from Eq. p.6p . one gets the simple differential equa- 
tion for the right electron asymmetry j/r = lO'^^gj^, 



dt 



(r + rB)j/R = Q, 



(5.1) 



where T — SFrl is the chirality flip rate 
Tb = 6(5'V47^2)2B2/100T^ and Q = 6 x lO'^ x 
g'^BYko/AOOiT^T^ come from the second (helicity) term 
in Eq. (j3.6p and from the first (diffusion) term in the 
same Eq. p.6p . The solution of Eq. ()5.ip obtained for 
strong and constant hypermagnetic fields Fb 3> F and 
By ~ const 



ynit) 



Q 



F + Fb 



-(r+rB)(t-to) 



F + Fe 



(5.2) 



gives the asymptotic growth of yn{t) up to yn{tEw) [here 
for its initial value ^eR(io) = 0], 



?/r(<ew) 

1 R 



10* 



9' 



') {tew) 



Tb 

: 0.32. 



(5.3) 



Here we put Fb^ew ^ 1 for strong fields, as well as 
substituted g'^ = e^/ cos^ 9w = 0.12 and ka/T^w = 10""^ 
for the case of Bq — 25.6. 

The bigger the wave number or the stronger the 
hypermagnetic field By — koY{t), the weaker a sphaleron 
influence on lepton (baryon) asymmetry is. In a wide 
region of wave numbers, the left lepton asymmetry ?/l 
remains negligible compared with the right electron y^ , 
Ij/lI J/R- This is due to the initial conditions p. lip 



in our scenario resulting in the right electron asymmetry 
?/Fi and the chiral anomaly 2/r — close to what was 
calculated in Ref. [l^. However, for a long CS wave 
ko <C fcmax, a small positive value j/l > evolving down 
to Tew [cf., Fig. EJd)] is sufficient to allow sphaleron 
transitions to wash out BAU. 



VI. DISCUSSION 

In the present work we have studied how the chiral 
asymmetry A/i = fie„ — fiei, ^ arises before EWPT 
in the scenario where, firstly, the initial right electron 
asymmetry S,en{to) = Men /To ~ 10~^° provides the gen- 
eration of BAU. Secondly, if we discuss the temperatures 
T < To = Trl ~ 10 TeV, chirality flip reactions enter 
the equilibrium. Then, at t > to, the violation of lep- 
ton numbers leads to a nonzero left electron asymmetry 
^ei,{t) = Mcl/T" 7^ 0. Note that we consider the lepton 
numbers violation due to Abelian anomalies and because 
of the presence of the SU(2)w anomaly. The generated 
left electron asymmetry results in change of the primor- 
dial right electron asymmetry S.eu_{t) and influences the 
BAU evolution. 

For large scales of hypermagnetic held (for a smaller 
value fco)j sphaleron transitions are more efficient to erase 
BAU since the amplitude of a competitive mean hyper- 
charge field decreases [when amplitude i3y(i) = koY{t) 
goes down], and therefore due to Abelian anomalies, an 
enhancement of the lepton number ceases. Of course, the 
bigger the seed hypermagnetic field By \ the bigger the 
lepton asymmetries and the baryon one. 

While our choice of CS wave as the simplest hyper- 
magnetic field configuration significantly simplifies the 
analysis of the lepton (baryon) asymmetry evolution, we 
should comment on some disadvantages that are appro- 
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FIG. 3. The normalized chemical potentials ?/r,l and the chiral anomaly parameter j/r — j/l versus f/tEW for B^> = 10" G. 
(a) The normalized chemical potential for Bo = 25.6 (solid line) and Bq — 2 x 10~^ (dashed line). The dash-dotted line 
corresponds to the asymptotic value of yn. — 0.32 calculated analytically in Eq. (|5.3[l . (b) The chiral anomaly parameter 
J/R — yL for -Bo ~ 25.6 (solid line) and Bo = 2 x 10"'^ (dashed line). The dash-dotted line corresponds to the asymptotic value 
of J/R — vl equal to 0.34. (c) The normalized chemical potential t/l for Bo = 25.6. (d) The normalized chemical potential yL 
for Bo = 2 X 10"^ (solid line) and Bo = 2 x 10"^ (dashed line). 



priate to such hypercharge field. 

First, a scale of hypermagnetic field for the chosen CS 
configuration is rather small. In order to get the BAU 
close to the observable value -Bobs(iEw) ~ 10~^° we fit 
the CS wave number of the order of /cq — lO^^fcmax, 
where fcmax = IO^^Tew corresponds to the maxi- 
mum wave number for the hypermagnetic field surviving 
Ohmic losses at EWPT time. The macroscopic long- 
range hypercharge field has the scale A = k^^, which is 
much bigger than the mean distance between particles in 
plasma T~^, and, on the other hand, is much less than 
the horizon size A <^ Ih- Here the horizon size, e.g., at 
the EWPT time Ih = {Mpi/1M^)T~^ = 
is much bigger than the scales — 10^/Tew and 

ko^ ~ (10^° - 10")/rEw applied in our plots. Be- 
cause of the arbitrariness of the z-axis direction cho- 
sen for the CS wave, such macroscopic fields are rather 
small-scale (random) fields and there is no anisotropy 



of medium. Thus, in order to get a necessary scale 
on the onset of galactic magnetic fields we should rely 
on an inverse cascade evolving after the EWPT for the 
Maxwellian fields [IJ] that originated from the hyper- 
charge ones in our causal scenario. Note that such in- 
verse cascade needs a significant amount of magnetic he- 
licity in order to operate efficiently. For the CS wave, 
the helicity density hy ^ kg J dtY^{t) decreases faster 
with lowering of the wave number ko than the energy 
density pif^ = = k'^Y'^{t)/2. In other words, the 

lower ko, the further away the CS wave configuration 
is from the maximum helical field obeying the relation 

( B] 

kohy'^^ — 2py . This circumstance should be taken into 
account for a more realistic continuous hypercharge field 
spectrum when the conservation of the global helicity 
governs a spread of helicities over different scales. 

In addition, our choice of CS wave does not appear to 
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be not realistic for transition of hypercharge field to the 
Maxwellian one during EWPT. It was shown in Ref. [l^ 
that being provided by the helicity conservation only 
a helical field such as that given by the 3D configura- 
tion of can penetrate the boundary wall separating 
symmetric and broken phases during the EWPT time 
T ~ Tew- The CS wave does not penetrate such surface 
of a bubble of a new phase even for a strong hypermag- 
netic field amplitude, which, in turn, provides possibility 
of EWPT of the first order for the present bounds on 
Higgs masses [20| . 

The evolution of the corresponding hypermagnetic he- 
licity for an arbitrary configuration of hypermagnetic 
fields before EWPT Hy = J d^ii{Y ■ By) has been re- 
cently studied in Ref. neglecting hypermagnetic dif- 
fusion. The following magnetic helicity evolution in hot 
plasma at temperatures T <C Jrw was analyzed in the 
same approximation in Ref. [27| relying on the model 
proposed in Ref. for the magnetic helicity parameter 
a ^ Gp. The new mechanism for the a-helicity parame- 
ter suggested in Refs. [Tl','T2i] may improve such estimates 
for primordial Maxwellian fields. 

To resume, we estimated the chiral asymmetry A/i = 
Men ~ arising just at EWPT T ~ Tew using the sim- 
plest configuration of the hypercharge field — the Chern- 
Simons wave — and taking into account both Higgs de- 
cays and inverse Higgs decays and sphaleron transitions 
as well. The evolution of lepton and Higgs asymme- 
tries were studied at temperatures Tew < T < Trl. 
The baryon and lepton asymmetries crucially depend on 
the CS wavelength kg^. The doubts whether it is possi- 
ble to protect the baryon asymmetry of the Universe in 
the symmetric phase by temporarily storing BAU in the 
asymmetry of the cr species are dispelled in the case of 
strong hypermagnetic fields. The washing out of BAU 
by the sphaleron transitions due to the involvement of 
left particles at T < Trl through Higgs (inverse) decays 
is not dangerous in a wide region of CS wave numbers. 
A strong seed hypermagnetic field i?Y ^ is needed to sup- 
port all such issues while dynamo amplification turns out 
rather negligible for the CS wave. The amplification of 
hypermagnetic fields by mechanisms beyond the SM, as- 
suming, e.g., a new pseudoscalar field coupled to hyper- 
charge topological number density 28] , is not considered 
here. Within the framework of the SM model, the 3D 
configuration of a hypermagnetic field seems to be much 
more productive for problems under consideration that 
remain as the challenge for us in the future. 
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Appendix A: The origin of the chiral anomaly for a 
Maxwellian field and the CS anomaly for a 
hypermagnetic field 

The chiral anomaly parameter Afi = {^e^i ~ Mel) 7^ 
leads to an additional contribution to the current in the 
Maxwell equation 

" 'dt 



V X B = (JcondE - 



TT 



A^B, 



(Al) 



where the last pseudovector current j = (aomAM/7i')B is 
coming, e.g., from the energy balance under chirality flip 
for massless particles (29l . l30l | 



/d3x(j.E) = ^^/ 



d^x(E ■ B) 



(A2) 



Note that acm and Ucond were defined in Sec. ID 

Using the Bianchi identity 9tB = —V x E, one finds 
from Eq. (jAip . neglecting in the MHD approach the dis- 
placement current 9tE — 0, the modified Faraday equa- 
tion Oil 



OB 

'dt 



V X B 



1 



-V^B, 



(A3) 



which governs the evolution of the magnetic field after 
EWPT at temperatures 10 MeV < T < Tew- 

Let us comment on the energy balance Eq. (|A2[) . Note 
that in the rhs of Eq. (|A2p the Adler anomaly for the 
right-handed electrons (positrons) 



TT 



E-B) 



and for the left-handed ones 1311 



16n 



^T^.T^'^-- — (E-B) 



being combined for the pseudovector — = (4'7''7^5') 
in uniform medium as 



dt 



(nR - ul) 



2ar. 



(E-B) 



(A4) 



defines the rate of the chirality flip per unit time per vol- 
ume and the corresponding energy cost (A/i/2)d(?T.R — 
Tii^)/dt = (j - E) for such effect. Hence, multiplying by 
A/i/2 and integrating Eq. (jA4| over volume, one derives 
the energy balance (|A2[) . While separating E in both 
sides in the integrand of (|A2p , one obtains the pseudovec- 
tor current j in the Maxwell Eq. ()Al|) . 

The CS anomaly term in the SM Lagrangian for the 
hypercharge field also leads to the pseudovector con- 
tribution in the Maxwell equation [T3| (here at tempera- 
tures Tew < T < Trl): 



a' ( 

V X By =(TcondEy H /Ig 

TT V 



-il 

"47r' 



(A5) 
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where the last term has polarization origin in the pres- 
ence of hypermagnetic fields in primordial plasma Q. 
Notice that the Adler anomaly used in Eq. (|A2|) corre- 
sponds to the difference of mean densities of the right- 
and left-handed fcrmions that is pseudoscalar [32] ?ir — 
riL = (^r^r) - = (*+7^^') with ^-r^l = 

(1 ± 7'^)*/2, while the CS term in Eq. (|X5t is given by 
the mean spin {'i+a.^^'^) = = M/^b, where 

pseudovector M ^ By is the magnetization of medium 
in a hypermagnetic field By- 

Let us stress the similarity of both chiral magnetic ef- 
fects leading to anomalous terms in Eqs. (jAip and (jASp . 
They are caused by a polarization mechanism provided 
by the main Landau level contribution to the additional 
current in the Maxwell equation (compare in Ref. (s^)- 



Appendix B: Lepton kinetics with Higgs bosons in 
the absence of hypercharge fields 

In this appendix we briefly discuss the lepton kinet- 
ics in the presence of Higgs bosons without hypermag- 
netic fields and Abelian anomaly and omitting sphaleron 
terms. In particular, we explain in detail the derivation 
of Eqs. dSi]) and 

a. Right electrons The reactions contributing to 
the right electrons dynamics are (i) the inverse decays 
grBl — > ip'--^^ and brP^ — > tp^-^'> and (ii) the direct decays 
(^(0) _j. brCl and (p^~^ — e^v]:. We remind that when 
particles annihilate, a "minus" sign appears in front of 
the Frl term. Taking into account that n^{o) — n^(-), 
as well as the equivalence ng-^ = n^L , one gets a factor of 
2 in the kinetic equation accounting for the two channels: 



n,J/» + 2rD^. (Bl) 



Here the Bose distribution for the Higgs doublet (p^ — 
{(f'^^K is given by the chemical potential jiQ — tJ-+ — 
— = ^A*(5) or fiQ = fj,^ for ip, while for the antiparticle 

ifi^ — ((yS*^^-*, i^'*'-') one gets /i_ = /ij-p^ = —fiQ. Note that 

equilibrium, e.g., in the reaction eRCL would 
correspond to the right relation of chemical potentials 
neglecting hypermagnetic fields as given by Eq. (I2.ip , 

A'eR — Mel — ~M0- 

b. Right positrons The reactions contributing to 
the right positrons dynamics are (i) the inverse decays 
grBl — ^ 9?^°^ and eRi/g — > (p''^\ as well as (ii) the decays 
tp^^^ — T' brCl and ip^~^^ — >■ cri^^ . Analogously to the right 
electrons case, we take into account that n^(o) = n^(+) 
and = n^i,. Finally, we obtain the following kinetic 
equation accounting for the two channels: 



d_ 

dt 



(^) = 2rRL[-neR - nej/s + 2rD- 



Subtracting Eq. 
count that L„„ ; 



(B2) 



B2I) from Eq. (|B1[) and taking into ac- 



] / s and Le-^ = [r 



and in our previous work [l3| if we omit the Abelian 
anomaly [see the first line in Eq. (13.11) in Sec. IIII| . 



dt 



2rRL(ieL - Leji) + 2rD[n^(o) - n (o)]/s. (B3) 



In equilibrium, dL^^ /dt — accounting for Frl — 2Fd, 
we get the correct relation (|2.ip . 

c. Left electrons We should take into account (i) the 
inverse decay crCl — > f''^^ and (ii) the decay, — >■ 
grCl, which give 



dt 



(B4) 



d. Left positrons In this case (i) the inverse decay 
eLGR — >■ (p'--^^ and (ii) the decay i^'"^ — blBr give the 
following contributions: 



^ " rRL(-nej, - UgJ/s + Fr 



dt \ s / - g 

Subtracting Eq. (jB5|) from Eq. (jB4[) one gets 



(B5) 



dL, 



dt 



rRL(-^eR ~ 



(0) 



■n^(o)]/s. (B6) 



In the equilibrium dL^i^/dt = accounting for Frl = 
2Fd, we get the correct relation for chemical potentials 

Men - AieL + = 0. 

Let us derive the Higgs boson kinetic equations. In 
Eq. (|B4[) . the boson ip^^^^ decays into the pair e^e^- 
Hence, the boson vanishes, which increases the popula- 
tion of Bl and cr (in kinetics of left electron bl). In the 
kinetic equation for ip^^'^ itself, such term enters with the 
opposite sign as — Fd'^<^(o) (boson disappears). Analo- 
gously, the inverse decay term should have the different 
sign in boson kinetics: it becomes -l-FRL(?T.eR + '^eL)/s 
since the pair crBl annihilates into tp^^^ increasing the 
population of the neutral bosons tp^^\ Therefore, one 
obtains from Eq. (IB4p . 

^(^)=rRUn.R+n.J/. + F,(-!^). (B7) 

Analogously from Eq. (jB5l) . changing sign on the right- 
hand side, one obtains the kinetic equation for the Higgs 
antiparticle (^j'^^h 



d /n^{o) 



d /? 
dt V 



rRL( )/s + rD - 



(B8) 



Subtracting Eq. (jB8|) from Eq. (jB7p and accounting for 
Fd = Frl/2, we derive the kinetic equation for the Higgs 
boson asymmetry [see Eq. p.2p in Sec. Illlj : 



d_ 

dt 



— [(n^(o) - n^(o))/s] 



one gets the equation similar to that derived in Ref. [1( 



— Frl S Le, — Le„ — 



2s 



(B9) 



11 



e. Left neutrinos For left neutrinos we account for 
(i) the inverse decay eR.i'g — 7> tp^~^^ and (ii) the decay 
(y9(+' — )■ BRi^g , which give 

:t7 ( — ) = rRLhng, - n,L]/s + Td^^. (BIO) 
dr V s / e g 

/. Left antineutrinos In this case (i) the inverse de- 
cay BRpg — )■ (/s*^"^ and (ii) the decay Lp^~'^ — > e^iv]^ con- 
tribute to the kinetic equation as 

d /?^pL 



one gets the well-known issue 9^Jr = that guaran- 
tees the conservation of the right electron current in the 
absence of hypermagnetic fields. On the other hand, sub- 
tracting Eq. (IC1|) from Eq. (|C2[) we get the violation of 
the left lepton current in the same non-Abelian fields, 



9 pi^f p 



(C3) 



Notice that here neutrino and electron currents are equiv- 
i^L = Jl/2, as seen from the following repre- 



alent j^^ 



_ (!!:EL) = TRLf-nej, - npLl/s + Td^^^. (Bll) scntation for the left field doublet L^ = (i>e , gl), 
at \ s / " s 



Subtracting Eq. (jBlip from Eq. (jBlOp . one gets 



dL„ 



dt 



Lei,) + ^i)[n^i.+) - n i-)]/s, (B12) 



where we took into account that L^^l — o^' "f^v^ — '^cl ■ 
Of course, = n^{o) with the chemical potential 

in Bose distribution for the doublet ip^ — ((^(+\ (^(°^), 
and n^(-) — n^[o) with the chemical potential — /iq for 
the c.c. doublet = ((p(-), (^(o)). 

For charged Higgs, which are described by Eqs. (|B10p - 
(|B12p , using arguments like in the derivation of Eq. (jB9p 
we obtain from Eqs. (jBlOl) and (IBlip the kinetic equation 
which is identical to Eq. (jB9l) . since n^p{+) = n^(o) and 



d_ 

At 



[(n^(+) - n^(-))/s] 



2s 



(B13) 



L = 



Va2 + 62 V 



where we may put a = — 6 = 1 for the isospin col- 
umn. Then, using the standard field operator in the 
Schrodinger representation 

+ <i!;(p)<(p)e-*"i. 

one finds the current asymmetry (x, t) = 

Tr[/5(t)?/;L7^V'L] = Ji^ "-^II' ^'^^'^^ P^^) ^-'^^ nonequilib- 
rium statistical operator obeying the Liouville equation. 
Here the currents 



d^p p'^ 



/('-'-)(p,x,0, 



Appendix C: SU(2)w anomaly and left fermion 
number violation 



are given by the Wigner distribution functions 



Let us use Eq. (12-174a) in Ref. '34] written for the 
pseudovector current of the one generation of massless 
fermions, = - j'^ = ■ipj'^j^ip, 



f 



(p,x,t) =Tr 



^ •'p-|-k/2,r';p-k/2, 

. k 



WR Jhi ^ 



(CI) which, in turn, are given by the distribution func- 
tions in the momentum representation f'^^r'-vr^^) ^ 



Here = '07''(1 -h7^)VV2 and = -f^)'ip /2 are 

the right and left fermion currents correspondingly [35| . 
Adding the equality (ICip with the anomaly for leptons 
of the first generation given by Eq. (11.12) in Ref. [lB| 
(see also [36.]). 



167r 



^i^r^^M-, (C2) 



Tr 
Tr 



p{t)blip)bAp') 

p{t)dUp)dr,{p') 



p'r' ;pr v 

for particles and fp^J.prW 
for antiparticles. 



The violation of the left lepton number L^^ in non- 
Abelian fields due to the SU(2)w anomaly (|C3p proceeds 
with the sphaleron transition probability Tgph as we used 
in kinetic equations (13.11) . 
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